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INTRODUCTION 

This short essay has two purposes: (I) A rapid review of classical and stochastic quan- 
tum Liouville equation and proof of that and some use of Liouville equation in statistical 
mechanics and (II) to present a short solution for a simple solvable model. The Liouville 
equation is valid for both equilibrium and non-equilibrium systems.lt is imperative to the 
proof of fluctuation theorem from which the second law of thermodynamics can be derived 
and also it is the key component of derivation of Green Kubo relation [ljfor linear transport 
coefficients such as shear viscositythermal conductivity or electrical conductivity. 
The dynamics of such composite systems is governed by a quantum-classical Liouville equa- 
tion for either the density matrix or dynamical variables which are operators in the Hilbert 
space of the quantum subsystem and functions of classical phase space variables of the 
classical enviroment.For more applications in quantum mechanics refer [2j. 



1 Corresponding author 
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I.DESCRIPTION OF CLASSICAL LIOUVILLE EQUATION AND ITS PROOF 

The Liouville equation is a cardinal equation of statistical mechanics. This equation de- 
picts the development of phase space distribution function for the conservative Hamiltonian 
system and also supplies a complete description of the system both at equilibrium and away 
from equilibrium. On the other hand Liouville equation is a persistent for the flux and 
meditate the phase space of micro-canonical ensemble of a N-particle system(NVE). Let 
p(q( N \p( N \ t) refer the phase space density. If we contemplate a volume element dp^dq^ 
,then p(q( N \ p( N \ t)dp( N ^ dq( N ^ acquire the number of distinctive points in volume element. 

It is a typical form of equation of motion of p(q( N \p( N \t). 

The number of system in the ensemble preserved the abidance equation for development 
of phase space density assumed as: 



dt 

dp „ 

implementing chain rule one obtains 

dp ^ , dp . dp / V"* ( dq ' ia l d Pic\\ /r)\ 

Ol 4=1,0=1,3 °9ta °t>xcx i=l, a =l,3 ° qta ° Pta 

In a Hamiltonian conservative system, the energy is kept up as a function of time and the 
time derivatives are measured by Hamiltonian equations of motion of classical mechanics. 
They are assumed by a set of coupled first order partial differential equations , 

*" = dp7 a (3) 

9H (AS 

dq ia 

i = l,2,...,N 
a = 1,2,3 

where qi a is the a — th component of the situation, qi of the % — th particle. 
If we use Hamilton's equation of motion we acquire: 

dqla = d 2 H ^ 

dq ia dq ia dp ia 
dp'ia _ d 2 H 
dP ia dp ia dq ia 
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thus the last term in equation (4) is indistinguishably zero leaving us with 

dp dp ■ dp ■ i 



up up . 

^ ^ dq ia qia dpi 



-p 



dp = dp dff dp dH 

dt 4^ dq ia dp ia dp ia dq ia 

The last equation is the classical Liouville equation. Some of scientists can use the classical 
Liouville equation symbolically : 

| = [*,d (8) 

Proof : consider an arbitrary " volume" u in the relevant region of the phase space and let 
the "surface" enclosing this volume be denoted by a. Then the rate at which the number 
of representative points in this volume increases with time is written as: 

|jU «» 

Where du = (d 3N qd 3N p) .On the other hand, the net rate at which the representative points 
"flow" out of u> (across the bounding surface a ) is given by: 

p(y • n)da (10) 

Here,f is the velocity vector of the representative points in the region of the (outward) unit 
vector normal to this element. By the divergence theorem(lO) can be written as: 



div(pv)du) (11) 
of course the operation of divergence here means the following: 

div(pv) = J2{ — (pq t ) + ^-(PP,)} (12) 

In view of the fact that are no " surfaces" or" sinks" in the phase space and hence the total 
number of representative points remains conserved 1 , we have, by (9) and (11), 

— / pdu = — [ div(pv)du (13) 



That is: 



J{^ + div(pv)}cLj = (14) 



1 This means that in the ensemble under consideration neither have any new members being admitted nor 
have any old ones being expelled. 
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Now, the necessary and sufficient condition that integral (14) vanishes for all arbitrary 
volumes u is that the integrand itself vanishes everywhere in the relevant region of the 
phase space. Thus we must have: 

Y t + d ™(pv) = (15) 

Which is theequationof continuity for the swarm of the representative points. Combining 
(12)and(15),we obtain: 



Q 3N r, 3N n . Q . 



The last group of terms vanishes identically because, by the equation of motion, we have for 
all i, 

dqt_ = d 2 H( y q i ,p i ) = d 2 H( y q i ,p i ) = _dpi_ 
dqi dqidpi dpidq { dpi 

Further, since p = p(qi,Pi, t),the remaining terms in (16) may be combined to form the total 
time derivative of p,with the result that: 

Equation (18) 2 embodies the so called Liouvillestheorem (1838). According to this theorem, 
the "local" density of the representative points, as viewed by an observer moving with a 
representative points,stays constant in time. 
Now we can understand : 

L = i[H,] (19) 
where i[H,] is a Liouville operator , symbolled by L and defined as: 

. r ^,dH d dH d s 

iL = \ ( ) (20) 

dpi dqi dqi dpi 



so we write: 



this has formal solution, 



I = < 21 > 



p(t) = exp(-iLt)p(0) (22) 



the classical Liouville equation is a highly non-trivial equation , where the momenta and 
coordinates of all the N-particles of system are in principle coupled with each other. 



We recall that the poissonbracket [p, H] stands for the sum: 

y^3W ldp_dH_ _ dp dH \ 
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II. PHYSICAL INTERPRETATION 

The quantity p(q( N \p( N \t)dp( NS) dq( N } is the probability that at a time t the physical 
system is in a microscopic state corresponded by a phase point lying in the trivial 6N- 
dimensional phase space element dp^N)dq^N). Hence the total number of systems in the en- 
semble is assumed by the integral over phase space of distribution , J p(p^ N \ q( N ~ ) )dp( N > dq^ . 

The Liouville equation is the 6iV-dimensional analogue of equation of persistence of an 
incompressible fluid. 

It explains the fact that phase point of ensemble are neither created nor destroyed. Thus 
we can interpret it as a conservation theorem similar to charge conservation role in Electro- 
dynamics. 



III. LIOUVILLE THEOREMS 



The initial value problem for the Liouville equation was first studied by Petrina and 
Gerasimenko [3| and then by Jiang j3],with its initial data in L 1 and L 2 , respectively. 
It has been shown by using the operator semigroup theory that this problem has a unique 
solution if its initial data belong to some subset of L 1 or L 2 . 



A:CAUCHY PROBLEM FOR THE LIOUVILLE EQUATION 

The Liouville equation is an evolution: 

j t f N (t) = -H»f N (t) (23) 

Where is the infinitesimal operator of the group Sj^(—t). Let us consider the liouville 
equation(23) as an abstract evolution equation in the Banach space L^(A) with initial 
conditions 

/at(O) = f N . 

THEOREM l.A. The Cauchy problem for the Liouville equation (23) has a unique solution 
in the space Ljy(A) of summable functions. It is given by the formula 

Mt) = S£(-t)&. (24) 

For initial conditions /jy G L° N (A) this solution is a strong one and for arbitrary G Ln(A) 
it is a generalized solution. 



PROOF. According to the well known result of functional analysis [103| the Cauchy problem 
for equation 

(5a (*)/iv)(^i, -,x N ) = 
hiX^t, x), X N (t, x)), (xi, x N ) G T N \W° N 

0, (gi, ...,q N ) G W N . 

has the unique solution in L^(A) which is given by (24). We shall demonstrate that in the 
formula (24) defines a generalized solution. For this purpose we consider the functional 

(tp,f N (t)) = / dx 1 ...dx N if{x 1 , ...,x N )f N {t,x 1 ,...,x N ) 
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Where tp is a continuously differentiable function with a compact support that vanishes like 
functions from L° N (A) in some neighborhood of forbidden configurations and dA. Since p 
is bounded and /jv(t) is summable ,the functional (<p, f N (t)) exists. Using (24) and the fact 
that S^(—t) is isometric we transform (tp, /tv(^)) to the form 



(f, /jv(*)) = (S A (t)ip, f N ) = j dxi...dx N (S A (t)ip)(x!, x N )f N (x u x N ). 

As in Theorem l.A we can show that the function (Sj^ (t)<p)(xi, xn) is differentiable with 
respect to t, and expression 



N 



as At — > uniformly with respect to (x%, £/v) belonging to compact sets. The functional 
(J2i=i Pi'^i^A 00^)> /&)> since tne function £\ =1 Pv&S^a (*)v)(^) is bounded. Therefore 
we can differentiate with respect to t under the integral sign in the functional (S^ (t)ip, f^).As 
a result we get 

, Nr. N o N o 

/*(*)) = (E p i ■ Qr(^(t)cp),f N ) = (s» (t) £ ^ ■ q-<p, ft) = (£ p * • ^ /*(*)) 

t=l * i=l ^ i=l 

and these equalities mean that the function = Sj^(— t)f^ is a generalized (weak) 

solution of the Cauchy problem for the Liouville equation (23). 

Theorem (l.A) holds also for a system of particles in dt v . We only have to replace (—t) 
and Hj!f by their counterparts S N (—t) and H N in L^. 

Let us construct the Liouville equa tion starting from the Hamiltonian dynamics of a system 
of elastic balls. As well known [1041 the class of pure states of Hamiltonian system of particles 



can be described by the distribution function 

N 

D NA {t, Xl ,...,x N )= Y[d{xik-X k {t,x )) (25) 

{ti, ...,%} k=l 

Where X(t, x°) is the phase point that is occupied by a system of N particle at time moment 
t if at the initial time moment t — the system was at phase point x° G §l vN x (A \Wn). 
Here the summation with respect to {ii,...,i N } is carried out over all the permutations of 
the indices i k £ (1, N). All possible states of the system considered can be described by 
the distribution function on the phase space T N (A): 

f N 
D NiA (t, Xl ,...,x N ) = / dxl...dx° N D N (xl,...,x N ) V Y[S(x ik -X k (t,x )) (26) 

Which is superposition of the pure states (25) of the system with a definite weight 
D^x®, x%). The last quantity is a known function of the distribution of particle co- 
ordinates and momenta that is specified at the initial time moment t = 0. Consequently an 
initial state of a system of iV elastic balls is , in general, described by the function 

f N 
D N)A (0,x 1 ,...,x N ) = / dx 1 ...dx° N D N (x° 1 , ...,x° N ) S~] FT 8(x ik - x° k ) (27) 



{h,...,i N } k=l 
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Which is equal to 

Av,a(0, xi, x N ) = N\D° N (xi, x N ). 

We note that the distribution function Djv,A(i) has the following properties: 

D NA (t, Xi,..., x N ) = D NA (t, Xi, qi,p*; q^p*; (28) 

\qi-qj\ = d,i^je (1,...,N), 
D NA (t,X!,...,x N ) = D NjA (t,xi, ...;qi,p*; ...,x N ),qi G dA, 
D NA (t,x 1 , ...,x N ) = 0, (q h ...,q N ) G W N . 

If we change the variables x° — > X(—t,x) in the integrand of (26), then after calculating 
the integrals we get the relation 

D Nt A(t, X!,..., x N ) = D NA (0, Xi{-t, x),..., X N (-t, x)). (29) 

From this formula and definition (2.3) of the evolution operator is follows that the state at 
an arbitrary time moment is connected with the initial state by the following formula: 

= (S^(-t)D NA (0))(x 1 ,...,x N ). (30) 

In order to construct an equation that determines the evolution of a state of a system of 
elastic balls we differentiate the right hand side of (30) with respect to time for D N a(0) £ 

From Theorem l.A it follows that 

d N 8 

—D NA (t,X!, ...,x N ) = -y^Pj. —D NA (t,x 1 , ...,x N ) (31) 

With the corresponding boundary condition in Poissonbracket and condition (28 ). The 
relation (31) is obtained as consequence of (29). From now on we shall treat (31) as an 
equation for an unknown function D NA (t,xi, ...,xn), which is the Liouville equation. We 
have to add to it initial and boundary condition (28), which are satisfied bjD^^t, x±, %). 



B:ON THE LIOUVILLE THEOREM 



Headword to study our problem by using the operator semigroup theoryequation (8) can 
be rewritten as 

Where A = -[H, .]. A can be in regarded to as either a distributional differential operator 
in L 2 (R) 2fc or a mapping of L 2 (R) 2fe into itself. 

THEOREM l.B .operator A is the infinitesimal generator of a contractive operator 
semigroup. 

In order to prove Theorem l.B , we'll first show the following lemmas. 



8 



LEMMA l.B . A is an increasingly-generated operator , i.e. , Re(Au,u) > for 
Vit G D(A) ,where (., .) denotes the internal product in L 2 (R 2N ). 
PROOF. By equation(8),we know that the for Vw, v G D(A), 

,a ^ [_^,dH dv dudH^ , 
(Au,u) = -J^-- r --) dq d P . 

By integrating partially and using the Hamiltonian equation ,the above integral can also be 
written as: 

f dH du du dH 
J jpl d Qi d Pi % dpi 

f — \^- ( 9H du du dH — - 

J ~i d 1i d Pi d 1i d Pi 

In particularly, 

(Au, u) = — (Au, u) 
Or equivalently, Re(Au,u) = 0, for Wu G D(u). 

LEMMA 2.B . A is a closed operator with dense domain in L 2 (R 2N ). 

LEMMA 3.B. XI + A is mapping of L 2 (R 2N ) on to it self , i.e., the rang R(XI + A) 
of the mapping XI + A is L 2 (R 2N ) , for any given A > , where / is identity mapping. 

For Vx G D(A). Then —A is the infinitesimal generator of S(t). 
Similarly, we can also show that 

THEOREM 2.B. operator —A is the infinitesimal generator of a contractive operator 
semigroup. 

According to the operator semigroup theory , by Theorem l.B , we define 
THEOREM 3.B. If u G D(A) then the initial value problem [(M + Au = 0), (u\ t=0 = u )} 
has a unique solution u = S(t)u G C\[0, +oo); L 2 (R™)) n C°([0, +oo); D(A)) 
where S(t) is defined by Eqs. [(S(t) = lim A ^ +00 S x (t)); (S x (t) = exp Akt ); (A x x = 
[X 2 (XI + AY - 1) - XI] X ); for^x e L 2 (R 2N )and, asX -> +oo; (A xX -> -Ax)] 

THEOREM 4.B. If w G .D(A fc ) , k is a positive integer , then the solution u of Eqs. 
[(^ + Am = 0), («| t=0 = u )} belongs to nJ =0 C fc - J '([0, +oo); D(A j )where, D(A°) = L 2 {R 2N ). 
For study their proofs refer 

The total time derivative of the phase space density is denoted by^. Then: 

* dt d qi dt dpi dV 1 J 

Where = , thus we provide the Liouville Theorem in classical Liouville equation : In 
conservative system the distribution function is constant along any trajectory in phase space. 
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The Petrina and Gerasimenko survey contains more and sound physical description 
in related to the operational approach to the solution of Liouville equation and non-linear 
ones(Bogolyubov). The Jiang argument constructed a weak group based on the non- 
differential operator version of the Liouville equation and the result as we can guess from 
the elementary group theory is a semi group represents a formal solution. AS Petrina and 
Gerasimenko discussed, this formal solution can be rewritten as a familiar Dyson series in 
terms of the normal generators of a complete group. However,since L 1 (R, 2K ) and L 2 (R 2k ) 
have not any mutual inclusion relation, the two existence results are different and the L 2 
existence theorem is not a theoretical scrutiny of the L 1 one. 



CrSOLUTION OF BOGOLYUBOV EQUATION 

The Bogolyubov equations for a system of elastic balls are the evolution equations 

jF{t) = —HF(t) + d 2 AF(t) (33) 

Where —H + d 2 A is the infinitesimal operator of the group U(t). Let us consider the Bo- 
golyubov equation (34) as abstract evolution equations in the Banach space L with the 
initial conditions 

F(0) = F° 



It can be show that |3| |104| . |103| that Cauchy problem (34) for the Bogolyubov equations 
has in the space L of sequence of summable functions a unique solution that is given by the 
formula 

F(t) = U(t)F° = exp( J dx)S(-t)exp{ J dx)F°. (34) 
In components (35) can be written as follows: 

F s (t, X%, X s ) — 2^ TTi TTj / dx s +i...dx s + n S S+n fe (— t, X\, Xg+n-k) x Fg +n (x±, x s + n j)35) 

n=0 fe=0 ■\ n ~~ )• J 

Each term on the right hand side of (36) is well defined, since the integrand is defined 
almost everywhere outside M® +n , and the series converges in the metric of space L s for 
t G] — oo, +oo[. 

THEOREM l.C. The Cauchy problem for the Bogolybov equations (34) has the unique 

solution in L that is given by (35). For initial conditions F° e L° this solution is a strong 

solution and for arbitrary F° e L it is a generalized (weak) solution. 

For study its proof refer [3[]. 

We note that the evolution equation 

±<p{t) = H<p{t) + d 2 A*p(t) 
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Are called the dual Bogolyubov equations 104 |. |l05l 



Obviously Theorem l.C holds also for the Bogolyubov equations that describe a system of 
particles moving in A: 

j t F A (t) = -H A F A (t) + d 2 AF A (t) (36) 

Where —H + d 2 A is the infinitesimal operator of the group U A (t) m(U(t) = 
exp( f dx)S(— t) exp(— f dx)). The solution of Cauchy problem for equations (37) is defined 
by 

F A (t) = U A (t)F% = exp( / dx)S A (-t) exp(- / dx)F A ) . (37) 

J A J A 

We not that the rigorous derivation of the Bogolyubov equations (34) ((37)) for a system 
of elastic balls, as an evolution equation in the space L(L(A)) of sequences of summable 
functions with the infinitesimal operator —H + d 2 A of the group U(t). Another method 
of justifying Bogolyubov equations consists in constructing these equations starting from 
Liouville equations (23). 



IV. THE STOCHASTIC LIOUVILLE EQUATION 

The transport equation is the stochastic Liouville equation ( SLE ).[5-9] 
It was shown. [Toj] . that the SLE can be re-expressed as presenting a /orma/trapping problem 
in a manner analogous to that employed in the theory of mutual annihilation of exciton's 
developed by one of the present authors. (l"D| 

The exact solution is particulary appropriate as a starting point for the study of the scatter- 
ing function relevant to experiments involving probe particles such as neutrons. The particu- 
lar feature that their solution. (Io| .possesses is the ability to address the degree of transport 
coherence of the moving particle which produces the scattering line shape. 
Before we explain section V we want to deal with spins which are in random motion. 



A:THE DENSITY MATRIX p(Q) 

The simplest example of such a type of motion is that of spins which jump back and forth 
among a number of sites with different chemical or magnetic environment. They (l2| .label 
the sites with the numbers v = l,2,...,n and define the distribution vector p v (t) as the 
vector of the probabilities of finding spins in the sites v at the time t. Another example is 
the case of translational motion. 

The position of a spin is then characterized by the spatial coordinate r and the distribution 
function is p(f, t) . The case of rotatory motion of molecules is very common. The orientation 
of the molecules is then defined by the Eulerian angles, fl All these cases in one formalism 
and choose Q as notation of the random coordinate. Thus the function p(Q,t) will denote 
the probability density of finding spins in the environment characterized by a particular Q. 
The spins can satisfactorily be described in terms of elementary magnets m,which can differ 
in direction but not in magnitude \m\. 
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The state of the whole spin system is then given by the combined probability density 
P(fl, fh), which we take to be normalized: 

P(Q,rh)dndm = l (38) 

From this function we can extract the following quantities;the molecular distribution func- 
tion: 

p(fi) = / P(tt,m)dm (39) 
The distribution function of the magnetization: 

f(m) = J P(tt,m)dtt (40) 

The average magnetic moment: 

(m) — / / mP(Q,m)dQdrh (41) 



The macroscopic magnetization of the sample: 

M = N(m) (42) 
Where N is the total number of spins per unit volume; and finally the quantities: 

m») = /«p(w (43) 

and 

M(Q) = N{m(n)} (44) 
Which may be called the Q — dependent magnetization density. 

These last two quantities have the following physical meaning. Suppose that we could measure 
separately the total magnetization of the molecules which have their orientations between 



Q, and f2 + cZO. Their measurement [12] , would then yield M(fi)dfLor suppose that we could 
measure the total magnetization of the sample, but with an Q — dependent weight w(ft).We 
would then measure j w(Q)M(Q)dQ.We can define (fn(Q)) in a slightly different way by 
writing: 

P(n,m) =p(n)q(n,m) (45) 

Where q(Q,m)is the normalized probability density of finding the magnetic moment at the 
value m, provided that we know that the molecule is in Q. We then have: 

(«(n)> = p(n)/fl 4 (n,*)*a (46) 

This notation emphasizes that (m(Q)) is the average magnetic moment of the molecules 
with orientation Q, multiplied by the probability of finding the molecules in fLThus 
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(m(f2)) contains information on the internal magnetic state of the molecules and on the 
molecular distribution in Q — space.lt is important to appreciate this point when the theory 
is applied to systems where p(£i) is not a uniform function, as is the case with chemical 
exchange between unequally populated sites, or with partially oriented molecules. 
It is not necessary to say, in the case of jumps the parameterf2is replaced by an index v and 
the integrations over Q became summations. 

We want consider molecular spin systems which are not adequately represented by a single 
magnetic moment m. The state of the spins has then to be described quantum mechanically 
by the wave function if) or the density matrix p. As for the density matrix,we should distin- 
guish between p of each separate spin system whose matrix elements p nm are the products 
a n a* m of the coefficient in the expansion of if) and its ensemble average, which we denote by 
p _1 ' 2 .If the spin system consists of a single spin / = |, there is a one to one correspondence 
between m and p, and between (m) and p~ 3 .However,if the spin system has more than two 
levels, this comparison can not always be made. Nevertheless it is often convenient to visu- 
alize the behavior of complicated spin systems through a model of magnetization. They 



utilized[12(] this analogy and assume a combined probability density P(p, f2)|13|],from which 



we drive the distribution function of the density matrices: 

f( P ) = J p(n, P )dn (47) 

The function / and P are defined in the space of all the possible density matrices p. This 
space is restricted to ps with matrix elements which are products of a n and a* m taken from 
normalized sets of coefficients a n . These are the matrices which satisfy the conditions of so 
called pure states. Not that the space of fa was also restricted to vectors with a constant 
modulus \m\. 

The determination of the functions P(Q,p) and f(p) is very difficult. 

We come to the main subject of introduction. Somewhat inexactly, call the fl— dependentspin 
density matrix 

p(fi) = J P P(n,p)dp (48) 

Which is the analogy of (m(Q)).lt may be used to compute the ensemble average of the 
expectation value of any physical quantity Q(Q) Which also depends on Q, 



(Q) = J tr[pQ(n)]P(Q,p)dQdp 

= J tr[p(n)Q(n)]dn (49) 

In most experiment performed on the spin system we measure a property Q which is inde- 
pendent of Q. It is then sufficient to know po i n order to calculate: 

(Q)=trQp ;p = J p(Q)dp (50) 

By analogy with equation(47) we can also write: 

p(n) = p(Q) J pq(Q, p)dp = p(Q)p n (51) 

Which emphasizes that p(fl) is the product of local ensemble average pn and the probability 
of finding molecules in Q. 
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B:THE DENSITY MATRIXp(i) AND SLE 

Let us now assume that the interaction of the spins are describe by an Q — dependent 
Hamiltonian H(Q),i.e. a Hamiltonian which is different for spins belonging to molecules 
with different points. How, then, does P(Q,p,t) develop in time if the molecules are moving 
randomly in Q — space! 

One starts to look at one particular spin system in the ensemble. Since this system belongs 
to a molecule which moves rapidly from one Q to the other, a time-dependent Hamiltonian 
H(t) is observed by the spins. H(t) is now written as the sum of a constant H and a time- 
dependent local Hamiltonian H'(t) with vanishing time average. 

The density matrix of spin system in question changes according to this Hamiltonian: 

(j t )p(t) = llp(t),H }+ l -l P (t),H\t)} (52) 

This equation is first solved by following a perturbation treatment, and then the ensem- 
ble average is taken. This yields the master equation. An essential point in the averaging 
procedure is that the correlation p(t) and H(t) is neglected. This is only permissible if the 
correlation times r of the matrix elements of H'(t) are so short that: 

\H'(t)\r<^fi (53) 

It is important to stress the difference between the density matrices p{t) and p(Q).p(t) is a 
pure state 12| and p(Q) is an ensemble average. Equation(53) is a stochastic equation which 



defines the stochastic process pit) in terms of the stochastic process H(t). More precisely by 
writing: 

(j t )p(t)= l Ji [p(t),H + H'{m}} (54) 

We see that process p(t) actually depends on the stochastic process VL(t). In most prob- 
lems dealt whit in magnetic resonance theory some model is assumed for the description of 
stochastic process Q(t).ln nearly all cases this is a stationary Markof fian process.lt is then 
assumed that the probability density p(Q, t) satisfies the equation: 

(^)p(n,t) = T P (n,t) (55) 

Where T is a time-independent Markof fian operator, operating on functions of fi.More 
generally, Q(t) is the projection of a Markof fian process, i.e. Q should be supplemented whit 
additional variables to form a complete set of random variables which make a Markof fian 
process. In order to retain a simple notation we assume that Q it self is a Markof fian 
process. Equation(56) fits well in formalism of section(IV.A) but it can be less directly applied 
to further development (53). Thus we follow Kubo's development [3] and write formally for 
the rate equation of P(f2, p, t) 

^P(0, p, t) = {-~\p, H(Q)] + T}P(Q, p, t) (56) 

This can be regarded as a composite Markof fian process. It is a coarse-grained descrip- 
tion of complete Liouville equation of the density of states of combined system of lattice 
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and spinsjitilizing the stochastic property T. Thus we may call it a stochastic Liouville 
equation (l3-18l . 



We first multiply equation (57) by p and integrate over p. This yields,with equation (49), 

^)p(n,t) = ~[p(n,t),H(n)] + rp(n,t) (57) 

This is the equation which we shall refer to as the SLE. 

In this derivation of SLE an important approximation has been made. It is assumed that 
the molecules execute their random motions regardless of the state in which the spins 
find themselves. Thus the reaction of spin system to its surroundings is ignored. In other 
words we neglected the energy exchange between the lattice and the spins. However, as 
Kubostated\13{ .this is permissible for instance when the temperature of the bath is suf- 
ficiently high compared with the possible energy exchange. Many examples in NMR or 
Mossbauer effects belong to this category because the reaction to the molecular motion of 
the bath is extremely small. Thus equation (58) has a wide range of application in line 
shape problems in magnetic resonance [181]. 

In the case of jumps the operator T is a matrix and the SLE takes the form: 

(j t )pu=^[Pu,H v ]+J2^P, (58) 

Where T va are reciprocals of mean residence times. 



V. SLE AND METHOD OF SOLUTION 

The stochastic Liouville equation in its simplest from is given by: 

Pmn ^[H, p\rnn Ot(\ &mn) Pmn (59) 

and describes the time evolution of density matrix p of the moving particle in the represen- 
tation of site states m,n. The system in which the particle moves is a crystal, i.e., possesses 
translational periodicity. The intersite interaction is H. The last in (60) describes the ran- 
domizing process whereby the off-diagonal elements of p decay, a being the rate at which 
this process occurs. Alternatively, a may be looked upon as the average rate of scattering 
among the band states of the particle. In the limit of no scattering, (60) describes coherent 
motion whereas, in the opposite limit it describes hopping or incoherent motion. 
The indices m, n are vectors in the appropriate number of dimensions. A different form of 
the SLE is: 

Pmn = ~i[H, p] - a(l - 5 mn )p mn + 25 mn '^('JmrPrr - JrmPmm) (60) 

r 

and describes in addition to the processes included in (60), a transport channel wherein the 
particle hops from sites r to sites m at rates 2^ mr . 

The SLE(60) can be looked upon as presenting a formal trapping problem in a 2d- 
dimensional space where dis the dimensionality of the system under analysis. We should use 
rj to denote the homogeneoussolution of (60), i.e., its solution in the last term: 



m—m',n—n' 

m'n' 



(61) 
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The quantities i\) are the density- matrix propagators, i.e., the solution of (60) for a = for 
the initial conditions p mn (0) = S m0 5 n0 . If we recast (60) as: 

Pmn ~\~ C^Pmn = ^[H, p\mn ~\~ Q^mnPmn (62) 

We see that a produces two perturbations on (62) : that caused by (ap mn ) on the left side 
and that caused by aS mn p mm on the right side. 

In the context of the dynamics of a hypothetical walker whose unperturbed motion is given 

by (62) , these two terms represent, respectively, an overall decay akin to the radiative decay 

of a moving exciton,[3], and a trapping or annihilation process, [7-9, 11-19], which takes place 

only when m = n,i.e, in a special trap-influenced region in the m, n space. 

The first of the two terms introduces a simple multiplicative factor in to the solutions (62) 

, but the second requires an analysis through the defect technique. fl9l|. 

Equation (63) takes the form: 



PrnnK^-) Vn 



m—m,',n—m' (e + a)p m ' m '(e)) (63) 



in the Laplace domain where e is the Laplace variable and tildes denote Laplace transforms. 
The case m = n gives: 

p mm (e) = f] mm (e + a) + a^(-0 m _ m / jTO _ m /(e + a)p m / m /(e)) (64) 

m' 

which involves only diagonal elements of the density matrix in the representation of m, n. 
We solve (65) through the use of discrete Fourier transforms. 
These are defined through relations such as: 

P k = ^ Pmm exp(ikm) (65) 

m 

where k is generally a vector and km a dot produce. 
The result is: 

1 — aip lc {e + a) 
It is straightforward to substitute (67)m(64). 

One then finds that the solution of the SLE(60) is given by (62) with the replacement: 
7 / \ T / \ a v-^r exp(zfc(s — r)) ~ ~ 

^ m _ m / jn _ n /(e) -> tp m - m >,n-n>(t) + T7 > ^[ 777 -\Vm-r,n-r{^ + ®)1p s - m > ,s-n' (e + "067) 

— , 1 — aw (e + a) 



T,s,k 



This result is exact and explicit. It is explicit in the sense that once one knows the ijj's ,i.e., 
the propagators of (60) or (63) in the absent of a, one can write down the solutions of (60) 
by following the prescription of (68) and the right-hand side of (62) , for arbitrary initial 
conditions. 

The practical usefulness of (68) depends on the simplicity, or lack thereof, of the quadrature 
problem involved in the inversions of the transforms. 
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VI. CLASSICAL LIOUVILLE THEORY AND THE MICROSCOPIC 
INTERPRETATION OF BLACK HOLE ENTROPY 

It has been computed that the entropy of a Schwarzschild black hole could be derived 
by finding a classical central charge of the Virasoro algebra of a Liouville theory and using 
the Cardy formula. This is done by performing a dimensional reduction of the Einstein- 
Hilbert action with the Ansatz of spherical symmetry and writing the metric in conformally 
flat form. Near horizon approximation the field equation for the conformal factor decouples 
becoming a Liouville equation. This computation is independent from a specific quantum 
theory of gravity model [20]. In order to understand the role of 2-D conformal symmetry in 



counting the black hole microstates let us recall some standard features of 2-D conformal 
symmetry. By introducing complex coordinates the flat metric can be written in the form 
ds2 = dzdz The infinitesimal generators of such transformations Gn = zn+l.z close a Lie 
algebra [Gn,Gm] = (n - m)Gn+m. 

In 2D CFT the stress energy tensor has only two components ,i.e., one analytic and the 
other anti-analytic: 

T = T ZZ 

f = f zz (68) 
d z T = 

d- z T = (69) 
We can therefore expand for example T in Laurent series: 

OO j 

r« = Ep£- (to) 

— oo 

In the quantum case the L n become operators and relatively to the commutator form a 
Virasoro algebra: 

C 

[L n , L m ] = (n - m)L n+m + — (m 3 - m)6 m - n . (71) 

This is a central extension of the algebra ([Gn,Gm] = (n - m)Gn+m). The central extension 
in the commutator algebra arises because of the normal ordering of the creation operators. It 
is a standard result for quantum CFT in 2-D that the asymptotic density of states for given 



L is completely determined by the Virasoro algebra by means of the Cardy formula 21 



p(L ) = exp(2 7 r V /^). (72) 

The entropy can therefore be calculated by using the logarithm of the Cardy formula S = 
ln p. 

A central extension of the conformal algebra ([Gn,Gm] = (n - m)Gn+m) can already arise at 
classical level in the Poisson algebra of the charges, as for example in the Liouville theory jl^. 



VII.APROXIMATE SOLUTION OF THE CLASSICAL LIOUVILLE USING 
GAUSSION PHASE PAKET DYNAMICS 



Their approach [23], to classical Liouville equation draws on the formal equivalence of 
this equation to the time dependent Schrddinger equation |24|. Gaussian wave packets,first 
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popularized by Heller 25|, |26|,for integrating the time dependent Schrddinger equation, have 
also been applied to the calculation of specteroscopic correlation functions using the von 
Neumann ( quantum Lioville ) equation 27H30 . 



Mukamel and co-workers observed that the equation of motion for a Gaussian approximation 
to the density matrix can be adapted to classical mechanics by replacing the quantum 
Liouville operator and phase space density distribution j3l[. The classical density distribution 
of each particle is represented by a single Gaussian phase packet. The many-body density 
distribution is expressed in a Hartree approximation as a product of distributions for each 



particle [32|, [33 



Equation of motion for the density distribution are derived for both constant energy and 
costant temperature dynamics and exact results are obtained for the time evolution in 
the free particle and harmonic potentials. We demonstrate the general applicability of the 
Gaussian phase packet dynamics method through two application. The first involves the 
calculation of equilibrium thermodynamic averages for a Lennard- Jones cluster and fluid 
using Gaussian phase packets. 



VIII. QUANTUM-CLASSICAL LIOUVILLE DYNAMICS OF PROTON AND 
DEUTERON TRANSFER RATES IN A SOLVATED HYDROGEN-BONDED 

COMPLEX 



They studied [34| an intermolecular proton transfer reaction in a bulk polar solvent of the 
form AH — B ^ A~ — H + B. The model under study, which was constructed by Azzouz and 
Borgis(35|. describes a hydrogen-bonded phenol (^4) trimethylamine (B) complex dissolved 
in methyl chloride. The proton transfer rate constant and kinetic isotope effect (KIE) have 
been computed for this model using a wide variety of techniques [3^-43]. The specific forms 
of the interaction potentials, parameter values used, and the remaining details of the model 
can be found in[36|, |44( and 45|. In a previous work, they calculated the proton transfer rate 
constant for this model with the AB distance constrained at Rab = 2.7 A. 



ArQUANTUM-CLASSICAL RATE THEORY 



The rate constant calculations are based on an expression for the time dependent rate 
coefficient of the proton transfer reaction A =f± B 



*Ub(*) 



1 



EEp 



n 



eq 



6a, 



ifiB 

dX x Re[N% a (X,t)W2 a (X, -£-)] 



(73) 



which is written in terms of a partial Wigner representation of the bath degrees of free- 
dom and a representation of the protonic degrees of freedom in adiabatic states 46[ . Here 
the bath phase space variables (coordinates and momenta, respectively) are denoted by 
X = (R, P),n e 2 is the equilibrium density of species A and j3 = is the inverse tempera- 
ture. In this partial Wigner transform representation [471] .the Hamiltonian H of the system 
is Hw = |E _|_ h(R), where h(R) is the protonic Hamiltonian in the field of fixed bath par- 
ticles with mass M. The adiabatic basis states {|a;i?)}are the solutions of the eigenvalue 
problem h(R)\a; R) = E a (R)\a; R), where {E a } are the adiabatic energies. In this expres- 
sion for the rate coefficient, N^ a (X, t)is the time evolved species £> (product) operator, while 
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W^ a (X, ^)is the spectral density function that contains all information on the quantum 
equilibrium structure. The solvent polarization reaction coordinate 48|, |49 . 



IX. STATISTICAL MECHANICS OF NON-HAMILTONIAN SYSTEMS 



The dynamics of Hamiltonian systems is characterized by conservation of phase space 
volume under time evolution 5(| , and this conservation of phase volume is a cornerstone of 
conventional statistical mechanics 5l|, [52J. Invariance of phase space volume under Hamilto- 
nian time evolution is the content of Liouvilles theorem for divergenceless flows 50_|, |53|. At 
a deeper level, conservation of phase space volume is understood to be a consequence of the 
existence of an invariant symplectic form in the phase space of Hamiltonian systems, and ap - 
plication of geometric methods and concepts from the theory of differentiable manifolds 54 - 



58|is essential for a fundamental description of classical Hamiltonian systems[5; 
Non-Hamiltonian dynamics, characterized by nonzero phase space compressibility j 



711 is relevant when we consider the statistical mechanics of thermostatted s yste ms 




Such systems arise in the simulation of ensembles other than microcanonical [76| , and in the 



treatment of nonequilibrium steady states [72|, |73|, |75|, [77 



Various homogeneous thermostatting mechanisms have been introduced to remove heat sup- 
plied by nonequilibrium mechanical and thermal perturbations. Phase space volume is no 
longer necessarily conserved, and for nonequilibrium steady states the phase space proba- 
bility distribution is found to collapse onto a fractal set of dimensionality lower than in the 
equilibrium case (72l. 73 . 75 . 78 . 79]. This phenomenon indicates a lack of smoothness of the 

8oL 8l| . The dynamical 



invariant measure in phase space in nonequilibrium steady states [7f 
evolution of the phase space distribution function for Hamiltonian systems is described by 
the Liouville equation |5( 
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The Hamiltonian equation is often considered be a special 
case of a so-called generalized Liouville eq uation (henceforth GLE) appropriate for systems 



eq i 

with compressible dynamics |53l. I62H71I. 82|, although the equation for the time-evolution of 



the Jacobian determinant in a general compressible flow given in Liouvilles original paper [53 



is, in fact, equivalent to the GLE|66|, |68|, [70j. A number of authors [6 2M7ll| have treated the 



statistical mechanics of non-Hamiltonian systems in terms of the GLE, and all have derived 
the result that the rate of change of the Gibbs entropy for non-Hamiltonian systems is the 
ensemble average of the divergence of the dynamical vector field (phase space compress- 



ibility). Steeb|67j, UM, [82 



applied the theory of Lie derivatives and differential forms to 
derive the GLE for both time-independent and time-dependent vector fields. Some explicit 
solutions to the Liouville equation were given, and the existence of singular solutions for 
systems with limit cycles (attracting periodic orbits) was noted j67[. The important paper 
by Ramshaw[7l[] gave a covariant formulation of the Liouville equation and of the entropy. 
It was noted that invariant measures (volume elements) associated with zero entropy pro- 
duction rate in non-Hamiltonian systems must be smooth stationary solutions of the GLE. 
Some limitations of the description of nonequilibrium steady states in terms of the GLE 
were discussed by Holian et al[83j. 

Tuckerman et al[84] (hereafter, TEA) have recently applied geometric methods from the the- 



ory of differentiable manifolds [55|-[58| , in particular the concepts of Riemannian geometry 57 



to the classical statistical mechanics of non-Hamiltonian systems [85j-[87j. TEA have ar- 
gued that, through introduction of so-called metric factors, it is always possible to define a 
smooth invariant phase space measure in non-Hamiltonian systems, even for nonequilibrium 



stationary states 85j, |86|. Moreover, the Gibbs entropy of the associated phase space distri- 
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bution function is found to be constant in time, just as for Hamiltonian systems. TEA also 
claim that previous formulations of the GLE (for exam ple, pap ers 

HQ) 

are m some way 

incorrect, incomplete, or at least coordinatedependent |85l437| . These claims have proved 
controversial 84, 



re- 



88H95 

By definition, a coordinate-free formulation in the language of differential forms (US 
moves any quest ion [8 71] concerning the coordinate dependence of any results obtained. The 
apparatus of differential forms is the appropriate machinery for treating the transformations 
of variables and volume elements arising in the dynamics of both Hamiltonian and non- 
Hamiltonian systems. For any region of phase space, the fraction of the ensemble inside the 
region is obtained by integration of the density n-form over the region. The density n-form 
can be written as the product of a volume formj55|, 56] (comoving volume element) and a 
phase space distribution function. The GLE, which describes the evolution of the phase 
space distribution function, then follows from the transport equation [67| . 
To determine the fraction of the ensemble in a given region, we simply need to count en- 
semble members. As we do not need a volume form to count ensemble members inside a 
prescribed region of phase space, the density n-form p is defined without reference to any 
particular volume form. Any result expressed in terms of p alone is therefore manifestly 
independent of the choice of volume form on the phase space manifold. This covariance 
with respect to the choice of volume form is the essential advantage of a description of the 
ensemble density in terms of the n-form p . 

When considering the phase space structure and dynamics of Hamiltonian and non- Hamil- 
tonian systems, the notion of distance associated with the familiar properties of Riemannian 
manifolds 55] is irrelevant. Volume forms provide exactly the construct needed, namely, a 
definition of volume without distance. Moreover, the Lie derivative of the volume form 
provides a definition of divergence without metric connection. The existence of an invari- 
ant volume form is important for simulation of equilibrium properties via non-Hamiltonian 
dynamics; if the dynamics preserves a given volume form (invariant measure) then, pro- 
vided the system is ergodic, and that all relevant constraints are taken into account [871]. 
phase space averages with respect to the invariant measure can be evaluated by computing 
long-time averages over a single trajectory. For all the equilibrium non- Hamiltonian systems 
discussed to date, a smooth stationary invariant measure can be found [871] . Nevertheless, 
any system with net attracting or repelling periodic orbits cannot possess a smooth invariant 
measure [96j . so that, contrary to the assumption made in 86], smooth invariant measures 
for arbitrary non-Hamiltonian systems do not exist. True invariant measures, which are in 
general singular (not absolutely continuous with respect to the usual volume element), can 
be defined in terms of infinite-time averages of the density n-form p t (8l| . 
If the metric tensor is stationary, then compatibility of the Riemannian structure with the 
dynamics requires the associated metric factor ^fg to be a time-independent solution of 
the GLE. In this case, the metric factor defines an invariant (Riemannian) volume form 



in the usual fashion [55 
dependent [8 



97|, [93 



On the other hand, if the metric tensor is allowed to be time- 
then compatibility requires the associated metric factor to be a non- 
stationary solution of the Liouville equation. The entropy defined with respect to the as- 
sociated time-dependent Riemannian volume form is then constant [§7j]; for nonequilibrium 
steady states, the underlying metric factor will become ever more nearly singular (fractal) 
at long times, so that this result is only of formal significance 
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X:SM ALL-TIME PARAMETER METHOD IN LIOUVILLE EQUATION 

Considering the formal asymptotic solution of Cauchy problem [J r.ffq , p, t) 



Lf(q,P,t),f(q,p,t)\ t =o = fo(q,p)] by the method of the small time space in[100| and de- 
fined by mapping r = 1 — exp(— st),s > Gennady Rudykh , Alexander Sinitsyn and 
Eugene Dulov transformed the initial infinite time interval R + onto a small finite one 
J A { r . o < T < 1}. 

This techni que is quite general and was applied to Liouville equation by G. Rudykh and 
A. Sinitsyn |l01 | at the earlier 80's of the 20th century. 

The small-time transformation is also well known to the applied mathematicians because 
it suits for numerical integration over the semi-infinite intervals. The obvious benefit for a 
such kind of transform lies in power series expansion over a time scale. Using finite interval 
[0, 1] one can pay the attention to the series coefficients and their convergence properties. 
In our particular transformed Cauchy problem becomes: 

(l-T)-^f(q,p,T) = hf(q,p,T) (74) 
f(q,P,r)\ T =o = fo(q,p) 

A solution of the Cauchy problem (75) in the small-time space was studied in the form of 
asymptotic expansion: 

oo 

f(q,P,r) = J2h(q,P)-r k . (75) 

fc=0 



Substituting (76) in to (75) and equating the coefficients for r , one obtains 102 



/I 1 n f) 

fk(q,p) = + ^([H(q,p),f k ^(q,p)} - £ ^{Qt{q,p) • /*-i(g,p)}) (76) 



Hence 



ft— l 

f k (q,p) = \l[H(1 + -L)]f (q,p) (77) 

Srv _, ST 

r=l 



With 



k = 2,3, .. 



k 1 1 1 1 

f[(l + —L) = a k _ x + -a k _ 2 L + ... + —doL^ 1 (78) 

SK S S 

r=l 



Where a = j^j, a k -i = 1- 



APPENDIX 



In this short section we present a concise proof for the theorem 2. A 
2.A:The family Sff (t),t e] — oo, +oo [, constitutes a strongly continuous one -operator group 
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of isometric operators in Ljv(A) whose infinitesimal operator is given on L^-(A) by the 
Poisson bracket with the Hamiltonian of system of non-interacting particles with bondary 
conditions. 

PROOF. Group properties of the family Sj^ (t) follow from the group property of a phase 
trajectory X(t,x)- For € L%(A) we show that: 

lim \\S%(t + At)f N - S»{t)f N \\ = 0. (79) 

At— >o 

Indeed, for f N 6 L° N (A) since 5^ (t) is isometric we have 

\\S%(t + At)f N -S%(t)f N \\ = \\SZ(At)f N -f N \\ = f dx\f N (X(At, X ))-f N (x)\- 

Since for sufficiently small At and fw(X(At,x)) is non-zero outside some neighborhood of 
forbidden configurations and the boundary dA of the domain A, where there are no particle 
collisions, we have f^(X(At, x)) ~ ► /at(x) uniformly with respect to x as At — > O.Thus we 
can carry out the lim At — >■ under the integral, hence the validity of (79) is proved. Since 
L° N (A) is dense everywhere in L7v(A),from (79) and the boundedness of the group S^(t) 
there follows the strong continuty of the group (t) in Ljy(A). 
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